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Abstract 

Electrical resistivity evolution during precipitation in Al-Sc alloys is modeled using 
cluster dynamics. This mesoscopic modeling has already been shown to correctly 
predict the time evolution of the precipitate size distribution. In this work, we 
show that it leads too to resistivity predictions in quantitative agreement with 
experimental data. We only assume that all clusters contribute to the resistivity 
and that each cluster contribution is proportional to its area. One interesting result 
is that the resistivity excess observed during coarsening mainly arises from large 
clusters and not really from the solid solution. As a consequence, one cannot assume 
that resistivity asymptotic behavior obeys a simple power law as predicted by LSW 
theory for the solid solution supersaturation. This forbids any derivation of the 
precipitate interface free energy or of the solute diffusion coefficient from resistivity 
experimental data in a phase-separating system like Al-Sc supersaturated alloys. 
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1 Introduction 



Resistivity measurements are usually a convenient way to follow precipitation 
kinetics in phase-separating systems. Indeed, assuming that resistivity is pro- 
portional to the solid solution solute content as it is usually done [1] , one gets 
direct access to the solute supersaturation. The so-called LSW theory (Lif- 
shitz and Slyozov [2] and Wagner [3]) as extended by Ardell [4] allows then 



emmanuel . clouet @cea. fr 



Preprint submitted to Elsevier 6 February 2008 



to deduce from these measurements key parameters like the solubihty, the 
precipitate interface free energy and the solute diffusion coefficient. 

Numerous such measurements exist in Al-Sc alloys [5-10]. They have been 
intended to characterize the precipitation of the AI3SC LI2 structure in alu- 
minum alloys. In a previous study [11], we used cluster dynamics to model pre- 
cipitation in Al-Zr and Al-Sc alloys. This mesoscopic modeling was shown to 
provide predictions of the precipitate size distributions in quantitative agree- 
ment with available experimental data [12-14]. It is therefore worth seeing 
if such an agreement can be obtained also with resistivity measurements. 
The purpose is mainly to see how predictive cluster dynamics can be, i.e. 
if one can obtain reliable information concerning the precipitates as well as 
the solid solution in the whole time range of precipitation kinetics and not 
only in the asymptotic limit of the coarsening stage like LSW theory does. 
Confronting cluster dynamics predictions with resistivity measurements will 
validate too the multiscale approach developed for precipitation in Al-Zr-Sc 
alloys [11,15-18]. The use of resistivity experiments to validate a multiscale 
kinetic modeling has already proved its efficiency for irradiated iron [19]. 

The first part of this article explains how resistivity can be modeled within the 
framework of cluster dynamics. A thorough comparison between the obtained 
model predictions and available experimental data is then performed. This 
allows us to conclude on the validity of the interface free energies which are 
obtained through such measurements. 



2 Cluster gas model of electrical resistivity 

Cluster dynamics models the phase-separating alloy as a gas of solute clusters 
which exchange solute atoms by single atom diffusion. Clusters are assumed 
to be spherical and are described by a single parameter, the number ngc of 
solute atoms they contain. The time evolution of the cluster size distribution is 
governed by a master equation which input parameters are the solute diffusion 
coefficient and the precipitate interface free energy. For Al-Sc alloys, these 
needed parameters were directly deduced from an atomic diffusion model [15, 
16,18] previously built for Al-Zr-Sc alloys. It should be stressed that, although 
none of these input parameters were intentionally fitted, cluster dynamics 
managed to reproduce reasonably experimental data on the time evolution of 
the precipitate size distribution. The reader is referred to Ref. 11 for a full 
description of cluster dynamics modeling and its application to Al-Sc alloys. 
We will only describe in the following how this technique can be adapted to 
model electrical resistivity of a phase-separating Al-Sc alloy. 
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2.1 Cluster contribution to electrical resistivity 



So as to use cluster dynamics to simulate the evolution of the electrical re- 
sistivity during the annealing of a supersaturated solid solution, we assume 
that the total resistivity of the cluster gas is given by the sum of each cluster 
resistivity. Doing so, we neglect any interference between clusters. Therefore, 
the input parameters of the modeling are the contributions of each cluster, 
i.e. pns^ for a cluster containing n^c solute atoms. Some calculations of this 
contribution exist in the literature (for a review, see Ref. 1). Most of them 
have been aimed to explain the increase of resistivity at the beginning of the 
precipitation kinetics in alloys hke Al(Zn), Al(Ag), Al(Cu) or Al(Cr), showing 
that this resistivity increase is associated with the apparition of small clusters 
(~ 10-20 A). 

All these calculations consider the elastic scattering of valence electrons by the 
perturbing potential due to the presence of the solute atoms composing the 
cluster. They show that solute atoms which are agglomerated in a small cluster 
can have a higher resistivity than the same isolated solute atoms because of 
Bragg scattering. Indeed, the smallest the cluster, the most relaxed the Bragg 
conditions for the electron scattering by the perturbing potential. This can lead 
to a high resistivity for small clusters. The scattering anisotropy increases with 
the cluster size and, as the cluster size gets of the order of the electron mean 
free path, the scattering becomes Bragg-like leading to a small contribution 
of the cluster to the resistivity. 

In Al-Sc alloys, such an increase of the resistivity in the beginning of the 
precipitation kinetics has never been observed experimentally [5-10]. Actually, 
it does not seem to exist in alloys where the solute is a transition element and 
the solvent a free electron like metal [20]. Therefore, precise calculations of 
each cluster contribution to the electrical resistivity may not be required for 
Al-Sc alloys. Moreover, these calculations are limited to small clusters as one 
assumes that the perturbing potential of the cluster does not modify the Fermi 
surface of the solvent and that electrons are scattered only once by the cluster 
(Born approximation). Such assumptions are reasonable only as long as the 
cluster does not contain more than a few solute atoms (~10) but do not really 
apply for larger clusters. On the other hand, for clusters with a stoichiometric 
composition, one can consider that the conductivity is infinite inside the cluster 
and that electrons are only scattered by the interface between the cluster and 
the matrix. For a sharp interface, this leads to a cluster contribution to the 
resistivity proportional to its cross section [1] and therefore to its interface 
area as clusters are assumed spherical. This leads then to pn^^ = Pi^sc^^^, and 
the electrical resistivity of the phase separating system at time t is simply 
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given by 

Pit)=PM + Pi E CnsMnsS' (1) 

nsc=l 

where C„g^(t) is the atomic fraction of clusters containing risc solute atoms. 
The time evolution of this cluster size distribution is governed by the cluster 
dynamics master equations (Eq. 2 and 3 in Ref. 11). In Eq. [T|, we have assumed 
that only clusters smaller than a critical size ng^. contribute to the electrical 
resistivity. We will see below which value this critical size has to be given in 
order to fit experimental data, p^i is the electrical resistivity of the solvent (Al 
in this case). It is temperature dependent and is known experimentally. The 
only remaining unknown parameter in Eq. [T] is the contribution pi of a solute 
monomer to the electrical resistivity. 

The increase of resistivity with Sc content has been measured at 77 K by 
Fujikawa et al. [6, 21] who give Spsc = 3400 nf2m per Sc atomic fraction. 
These measurements have been performed in under-saturated and therefore 
very dilute Al-Sc solid solutions. One can reasonably assumes that most of the 
solute are monomers and ignore larger cluster|3, with the consequence that 
pi = ^Psc = 3400 nf2m. Assuming that Matthienssen rule [1] is obeyed, this 
quantity does not depend on the temperature. Different measured values of 
6psc ranging from 3000 to 8200 nQm per Sc atomic fraction can be found in 
the literature (for a review, see Tab. 2 in Ref. 10). Nevertheless we find that 
the value measured by Fujikawa et al. [6,21] is the one that leads to the best 
agreement between the simulated and the experimental resistivity variations 
during precipitation kinetics. 

2.2 Critical size 

As already recalled above, the input parameters of cluster dynamics, i.e. the 
precipitate interface free energy and the Sc diffusion coefficients, were directly 
deduced from an atomic diffusion model [11]. The only left parameter needed 
to calculate electrical resistivity is the critical size n^^ appearing in Eq. [TJ So 
as to determine this boundary between clusters which are contributing or not 
to the electrical resistivity, we compare our predictions using different values 
for rig^ with experimental data (Fig. [1]). If one assumes that only monomers 
are contributing to the electrical resistivity (rig^ = 1), cluster dynamics do 
not reproduce experimental data. Indeed the simulated resistivity decreases 
too quickly during the precipitation kinetics. So as to obtain the plateau ex- 
perimentally observed at the beginning, one has to take into account the con- 

^ The maximal solubility limit of Sc in aluminum is 0.288 at.%. We have checked 
with our cluster gas model of electrical resistivity in Al-Sc that only the monomer 
contribution is relevant for such a low nominal concentration. 
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tributions of small clusters (nsc < 10), especially the dimer one. With these 
contributions, the simulated resistivity decreases at the right time, but the 
final plateau experimentally observed is not reproduced. Indeed, cluster dy- 
namics predicts that resistivity decreases to reach a value close to the one 
corresponding to the equilibrium solid solution, whereas experimental data 
show an excess resistivity. The final plateau corresponding to this excess resis- 
tivity can only be obtained if one considers that all clusters are contributing 
to the resistivity (ng^. = ocO])- With this value of the critical size, the exper- 
imental time evolution of the resistivity is reasonably reproduced by cluster 
dynamics. 

So as to better understand the evolution of the electrical resistivity we mon- 
itor the cluster size distributioilEl and each cluster contribution to resistivity 
(Fig. [2]) at different times of the precipitation kinetics for the supersaturation 
and the annealing temperature corresponding to Fig. [TJ At the beginning of 

^ Cluster dynamics equations can be solved only for a finite number of classes. 
The infinite size in the simulations corresponds to a maximal size which can be as 
large as ~ 10^^ solute atoms. We check that the concentration of clusters having 
this maximal size does not evolve during the simulation. If this is not the case, the 
maximal size is increased. 

^ Looking at the cluster size distribution, one should notice the abnormally high 
concentration of clusters containing 8 solute atoms. This is due to the low interface 
free energy associated with this size because of the existence of a compact cluster 
corresponding to a cube for this size [15]. 
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Fig. 1. Time evolution of the resistivity p experimentally observed [10] and deduced 
from cluster dynamics simulations for a solid solution of composition x^^ = 0.12 at.% 
annealed at 230° C. Different critical size ng^ have been used in cluster dynamics 
simulations (Eq. [1]). 
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the kinetics, {t — 10^ or 10^ s), only small clusters (nsc < 10) are present. 
Therefore, one can considers that resistivity only arises from these clusters. 
At the very beginning (t < 10^ s), one can even only take into account mono- 
and dimers. As precipitation goes on, the solid solution becomes depleted and 
thus the atomic fractions of small clusters decrease. For t > 10^ s, one can 
neglect the contributions of these small clusters to resistivity. On the contrary, 
as the precipitating phase appears, more and more large clusters are present 
in the system and only these large clusters contribute to the resistivity. For 
instance, for t — 10^ s, almost all the resistivity is due to clusters containing 
between 10 and 50 solute atoms. Thus large clusters are responsible for the 
excess resistivity observed at the end of the precipitation kinetics. As we are 




1 10 100 1000 

"Sc 

Fig. 2. Time evolution of the cluster size distribution (C„g^), of the cluster con- 
tribution to electrical resistivity {nsc^^^CnsJ and of the cumulative contribution 
iJ2n=i Cnn'^^^) in a solid solution of composition = 0.12 at.% annealed at 
230°C. Resistivities have been normalized by pi. 
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in the coarsening stage by this time, the concentrations of these larges clusters 
evolves very slowly, leading to a nearly stable resistivity. 



3 Comparison with experimental data 



Table 1 

Summary of the experimental conditions (temperature, concentration x^^, supersat- 
uration x^^/x^) corresponding to the available data used for the comparison with 
cluster dynamics simulations of the electrical resistivity variations in Al-Sc alloys. 





T^PTTi Ti Pi* n "1" n TP 

J- dll^J KjL CI t Lll C 




Sm m p v t; a f 1 1 T" fi I" 1 o Ti 


V\ fA'\TQ.f^i' anH RiT'iiTn 111 




19 at % 


'^5'^n 






1 2 at % 


874 




27n°C 


1 2 at % 


2fi5 




330° C 


0.12 at.% 


60 


Zakharov [8] 


250°C 


0.24 at.% 


942 




300° C 


0.24 at.% 


243 




350°C 


0.24 at.% 


77.7 




400° C 


0.24 at.% 


29.5 


Watanabe et al. [9] 


400° C 


0.17 at.% 


20.9 




450° C 


0.17 at.% 


9.05 


Jo and Fujikawa [6] 


260°C 


0.09 at.% 


264 




260° C 


0.15 at.% 


440 




300° C 


0.09 at.% 


91 




300° C 


0.15 at.% 


152 




370° C 


0.09 at.% 


19.4 




370° C 
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We now compare cluster dynamics predictions with all the different experimen- 
tal resistivity data available in the literature. All the experimental conditions 
corresponding to these data are given in Tab. [H When doing this comparison, 
we use the critical size nsc = oo, i.e. we assume that all clusters contribute to 
the electrical resistivity. 
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3.1 Al - 0.12 at.% Sc 



R0yset and Ryum [10] studied an alloy with a nominal compositioill] Xg^. = 
0.12 at.%. They followed the resistivity evolution during the precipitation 
kinetics for different annealing temperatures between 190 and 470°C. All re- 
sistivity measurements were performed at room temperature for which the 
pure Al resistivity was measured to be p\i = 27.0 nQm. For temperatures 
lower than 330°C, they observed that precipitates remain coherent and that 
precipitation is homogeneous. One can therefore compare their resistivity ex- 
perimental data with cluster dynamics simulations (Fig. [3]). A good agreement 
is obtained, especially for the lowest temperatures. Indeed, simulations man- 
age to reproduce the resistivity fast decrease during the nucleation and growth 
stage as well as the slower variation which follows during coarsening. Never- 
theless, for the highest temperature (330° C) which corresponds to a lower 
supersaturation (Tab. [1]), simulations appear to be too fast compared to ex- 
perimental data. The shape of the curve is correct, but there is a timescale 
shift between the simulated and the experimental resistivity. 

R0yset and Ryum [10] used their resistivity measurements to follow the Sc 
transformed fraction. In their treatment, they estimated the equilibrium solid 
solution from extrapolation of resistivity measurements in the coarsening stage 

^ They actually measured by EDS analysis a slightly higher composition, Xg^ = 
0.138 at.%, but cluster dynamics simulations obtained with this measured concen- 
tration do not significantly differ from those obtained with the nominal composition. 
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to infinite time. They found tliat tliis method gave erroneous solvus estimates 
at low precipitation temperatures. Our simulations (Fig. [3]) show that this is 
because of the excess resistivity associated with large clusters in the coarsen- 
ing stage. For high annealing temperatures (330°C for instance), the measured 
resistivity is close to the equilibrium one, but for the lowest temperatures the 
difference cannot be neglected. When taking into account this excess resistiv- 
ity one gets Sc transformed fractions different from those obtained by R0yset 
and Ryum (cf. appendix This explains why these authors deduced a wrong 
solubility limit for Sc in aluminum (Fig. 11 in Ref. 10), which they pointed out 
themselves, from their resistivity measurements. As cluster dynamics param- 
eters were obtained so as to reproduce the experimental solubility limit [11] 
and as our simulations manage to reproduce R0yset and Ryum resistivity mea- 
surements, we can conclude that these measurements completely agree with 
the assessed Sc solubility limit [22]. 



3.2 Al - 0.24 at.% Sc 




10 ^ 10° 10^ lO"^ 10^ 

t (s) 



Fig. 4. Time evolution of the resistivity p experimentally observed [8] and de- 
duced from cluster dynamics simulations for a solid solution of composition 
xl^ = 0.24 at.%. 

Zakharov [8] measured the electrical resistivity variations in an alloy contain- 
ing 0.24 at.% Sc annealed at different temperatures between 250 and 400°C. 
He did not specify at which temperature he performed these measurements, 
nor the resistivity of pure Al. We assume in our simulations a resistivity 
= 28.4 nf2m, corresponding therefore to a higher measurement tempera- 
ture or a lower aluminum purity than R0yset and Ryum measurements. 
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Cluster dynamics reproduce quite well these experimental data (Fig. Hj). For 
the lowest annealing temperatures (250, 300 and 350°C), the agreement is 
perfect in all the different stages of the precipitation kinetics. For the highest 
temperature (400°C), the simulated resistivity decreases too fast compared to 
the experimental one in the nucleation - growth stage. One should notice that 
the predicted evolution for this temperature is really fast and that the resistiv- 
ity drop appears at a time {t ~ 10 s) which is too small to be precisely observed 
experimentally. Nevertheless, even for this temperature, simulations manage 
to reproduce the slow decreasing of electrical resistivity during coarsening. 



3.3 Al - 0.17 at.% Sc 
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Fig. 5. Time evolution of the resistivity p and of the mean precipitate radius r 
experimentally observed [9] and deduced from cluster dynamics simulations for a 
solid solution of composition Xg^ = 0.17 at.%. The cutoff radius used to define 
visible precipitates in cluster dynamics is r* ~ 0.75 nm. 

Watanabe et al. [9] studied an aluminum alloy with 0.17 at.% Sc where they 
followed the resistivity during precipitation and measured the mean precip- 
itate radius using transmission electron microscopy. All measurements were 
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performed in the coarsening stage, after the resistivity had dropped because 
of precipitate nucleation and growth. As a consequence, only a partial compar- 
ison can be made with cluster dynamics (Fig. The reasonable agreement 
obtained between the simulated and the experimental resistivities is not really 
significant as resistivity does not vary too much in the observed time range. 
Due to the high temperatures of study (T > 400°C), the experimental as well 
as the simulated resistivities during coarsening are already close to that of the 
equilibrium solid solution. 

Cluster dynamics manage to reproduce too the variations of the precipitate 
mean radius (Fig. For radii greater than ~ 15 nm, Watanabe et al. ob- 
served that semi-coherent precipitates coexist with coherent ones. It should be 
pointed out that coherency loss is not taken into account in our simulations 
which handle only coherent precipitates. Despite this assumption, cluster dy- 
namics predictions are not too bad. This may indicate that the fraction of 
semi-coherent precipitates is small or that these precipitates have an interface 
free energy not too far from the one of coherent precipitates. 

3.4 Al - 0.09 at.% Sc and Al - 0.15 at.% Sc 

Jo and Fujikawa [6] followed the resistivity variations in two different alu- 
minum solid solutions containing 0.09 and 0.15 at.% Sc. The simulated evo- 
lutions of the electrical resistivity do not reproduce well their measurements 
(Fig. [H]). The resistivity drop associated with the nucleation and growth stage 
appears to occur too fast in simulations. Moreover, the slow decrease of the 
excess resistivity during coarsening is not as well reproduced as with other ex- 
perimental data. Here, only a semi-quantitative agreement could be obtained. 

3. 5 Summary 

Some more experimental data on resistivity measurements in Al-Sc alloys are 
available in the literature, but they cannot be used for a comparison with our 
simulations. Indeed, Nakayama et al. [7] studied an alloy containing 0.138 at.% 
Sc aged at 250, 300 and 350° C. But they normalized their measurements and 
did not specify the measuring temperature nor the corresponding pure Al 
resistivity, thus forbidding any use of their data. Drits et al. [5] studied an 
alloy the Sc composition of which is larger than the Sc solubility limit and 
even the eutectic composition. Some primary AI3SC precipitates may have 
appeared during the solidification and the alloy could not have been homog- 
enized. Therefore, one can doubt that this alloy initial state corresponds to 
a homogeneous supersaturated solid solution as assumed by our simulations. 
Nevertheless, all the previously reviewed experimental data already allow to 
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Fig. 6. Time evolution of the resistivity p experimentally observed [6] and deduced 
from cluster dynamics simulations for two solid solutions of composition Xg^, = 0.09 
and 0.15 at.%. 

draw main conclusions about the ability of cluster dynamics to predict elec- 
trical resistivity evolution during precipitation kinetics in Al-Sc alloys. 

Except for Jo and Fujikawa's data [6] for which only a poor agreement could 
be obtained, cluster dynamics manage to reproduce reasonably well experi- 
mental measurements of electrical resistivity [8-10]. For the highest annealing 
temperatures corresponding to the lowest supersaturations, a time shift could 
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appear between the simulated and the experimental resistivity, the time evo- 
lution predicted by cluster dynamics being too fast. Nevertheless, in all cases, 
the global shape of the evolution is correctly predicted. One can guess that 
a better description of the resistivity drop during precipitate nucleation and 
growth may be obtained if one would calculate precisely the small cluster 
contributions. But our simple model assuming a cluster contribution to the 
resistivity proportional to its area already leads to good predictions. Most 
importantly, cluster dynamics manage to reproduce quantitatively the slow 
resistivity decrease during coarsening. This decrease arises from the evolu- 
tion of the large cluster distribution and the surface dependency assumption 
appears to be correct for these clusters. 

An improvement of the agreement between cluster dynamics simulations and 
experimental data may arise too from a more precise calculation of the emis- 
sion and condensation rate coefficients appearing in cluster dynamics master 
equation as suggested by Lepinoux [23,24]. But, doing so, we will lose one of 
the main advantages of our approach, i.e. the ease with which the parameters 
of the mesoscopic modeling can be deduced from a very limited number of 
input data. In view of its simplicity, the ability of our approach to reproduce 
resistivity measurements appears more than reasonable. 



4 Discussion 



Resistivity measurements are often combined with precipitate size determina- 
tion using transmission microscopy so as to deduce from experimental data 
alloy parameters like the solute diffusion coefficient, the solubility limit, and 
the precipitate interface free energy. This can be done with the help of the 
LSW theory [2,3]. Indeed, Lifshitz and Slyozov [2] and Wagner [3] show that 
the precipitate mean radius varies linearly with the power 1/3 of the time in 
the coarsening asymptotic limit. Ardell [4] then extended the theory to predict 
the variation of the solid solution concentration. Applying his results [25] , the 
scandium concentration of the aluminum solid solution at time t should be 
given by 

Xsc{t)=4l + {Kt)-'/\ (2) 

where is the scandium solubility in aluminum and the rate constant k i 

_ Dsc ( kT \\ 



^ When using results obtained by Calderon et al. [25] for non pure precipitates, we 
assume an ideal solid solution and we use the fact that Xg^ << Xg^. 
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Dsc is the scandium impurity diffusion coefficient in aluminum, Xg^ = 1/4 is 
the Sc atomic fraction in the precipitate, a is the infinite radius hmit of the 
interface free energy between AI3SC precipitates and aluminum, and Q = a^/A 
is the mean atomic volume corresponding to one lattice site (a = 4.032 A for 
Al). Recently, Ardell and Ozolins [26] showed that the solute concentration 
can vary as the inverse square root of the time, and not the cube root like in 
Eq. [21 in the case where the precipitates present a ragged interface. This is not 
the case in Al-Sc alloys as the AI3SC precipitate interfaces are rather sharp as 
shown by our atomic simulations [15,16]. Therefore, one expects Eq. |2]to hold 
for this system. 

To make use of resistivity measurements, one usually assumes that resistivity 
depends linearly on the solid solution concentration. With the help of Eq. [2] 
and [3l one can then get information on the desired parameters, i.e. Dsc, 
or a. We previously saw that this linear relation between electrical resistivity 
and solute concentration does not hold in a phase-separating supersaturated 
Al-Sc solid solution as the resistivity has to be proportional to the cluster 
mean section (Eq. [T]). It is worth looking at the error due to the fact that 
one identifies the solute concentration with the cluster mean section when 
exploiting resistivity measurements. 

To do so, we need first to define the solid solution concentration in the cluster 
dynamics simulations. This is not so easy as this modeling technique describes 
the phase-separating alloy as a gas of solute clusters. Therefore, it does not 
differentiate between the solid solution and the precipitates at variance with 
precipitation classical descriptions like LSW theory (cf. Ref. 27 for a better 
understanding of the differences between cluster dynamics and classical theo- 
ries) . Nevertheless, one can discriminate the solid solution and the precipitates 
with the help of a threshold size ng^. Below this size, clusters represents fluc- 
tuations in the solid solution and above it they represent stable precipitates. 
The solid solution concentration is thus given by 

"Sc 

XScit) = ^ScCnM (4) 

"Sc = l 

In a supersaturated solid solution, one natural choice for this threshold size 
is the critical size rig^.. In cluster dynamics, this is the size for which the 
condensation rate Pnsc is equal to the emission rate ans^- Below this size, Pnsc 
is smaller than a„3^ and clusters have more chance to re-dissolve themselves 
than to grow. This definition works as long as we do not enter in the coarsening 
stage. A minimum for a size ngj,™ then appears in the cluster size distribution 
(ngj,™ = 3 for t = 10^ and 10* s in Fig. [2]). Once the critical size ng^ gets higher 
than ng^™, the quantity of matter contained in clusters smaller than ng^ begins 
to increase artificially because of small precipitates which become unstable. 
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We then choose the following definition for the threshold size: 

• rigj = rig^, as long as the cluster size distribution does not show any mini- 
mum, 

• ^Sc ~ ('^sc '^Sc™) otherwise. 
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Fig. 7. Change as a function of in the Sc concentration xsc{t) (Eq. and in 
the normalized resistivity [p{t) — p\]) /Spsc (Eq. [T]) simulated by cluster dynamics 
in an aluminum solid solution containing 0.12 at.% Sc annealed at 270°C and 330°C. 
The linear relation predicted by LSW theory (Eq. [2]) is shown for comparison. 

Using this definition of the threshold size, we monitor the variations of the 
solid solution concentration (Eq. H]) in clusters dynamics simulations for the 
annealing at 270 and 330°C of a solid solution containing 0.12 at.% Sc (Fig. [7]). 
The asymptotic behavior of xsc clearly obeys a linear dependence with t^^^^ 
as predicted by LSW theory. Using the solid solubility Xg^, the diffusion coef- 
ficient Dsc and the interface free energy a which the simulations rely on (cf. 
Ref. 15 and 11 for a full description of the way clusters dynamics parameters 
were obtained), we can calculate coefficients entering Eq. [2l The comparison 
with cluster dynamics simulations (Fig. [7]) shows that LSW perfectly manages 
to predict the asymptotic behavior of the solid solution concentration Xsdt)- 
On the other hand. Fig. [7] shows that the asymptotic behavior of the nor- 
malized resistivity (pit) — p^^i) /^Psc differs from the one of the solid solution 
concentration. Both quantities tend to a limit close to the scandium solubility 
Xg^ in aluminum, but the resistivity cannot be assumed to vary linearly with 
If one does so, the proportionality coefficient will differ from the one 
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predicted by LSW theory (Eq. [3]). As a consequence, if one assumes that the 
resistivity is proportional to the sohd solution concentration and uses LSW 
theory to deduce alloy properties from experimental data, the obtained solid 
solubility will be correct, but not the proportionality coefficient k. In particu- 
lar, one cannot deduce from resistivity measurements any reliable value for the 
precipitate interface free energy a or the solute diffusion coefficient D^c- The 
lower the annealing temperature, the bigger the error on these parameters. For 
instance, fitting with Eq. [2] and [3] the normalized resistivity evolution in Fig.[71 
one would get a = 2200 instead of 113 mJ.m~^ at T = 330°C. This clearly 
illustrates that resistivity measurements cannot be used to determine interface 
free energy in a phase-separating alloy as long as precipitates contribute to 
resistivity. 

Of course, this conclusion does not hold anymore when this precipitate contri- 
bution cancels. Indeed Watanabe et al. [9] manages to observe a linear varia- 
tion of the resistivity with This behavior was obtained once the precipi- 
tates became incoherent. Therefore, one can reasonably assume that incoher- 
ent AI3SC precipitates do not contribute anymore to resistivity. The interface 
free energies they deduced from their resistivity measurements (230 mJ.m~^ 
between 400 and 450°C) is higher than the one used in our simulations (119 > 
a > 105 mJ.m~^ between 200 and 500°C) which was deduced from the atomic 
diffusion model of Ref. 15, 18. This is in agreement with the fact that inco- 
herent precipitates should have a higher interface free energy than coherent 
ones. 



5 Conclusions 



Cluster dynamics has been shown to be able to reproduce electrical resistivity 
variations during precipitation in Al-Sc alloys. One should recall that the only 
input parameters required by this modeling technique are the solute diffusion 
coefficients, the precipitate/aluminum interface free energy and the resistivity 
increase per solute atom. None of these parameters was adjusted to fit the 
experimental data. Indeed, the experimental data measured by Fujikawa et 
al. [6, 21] was used for the resistivity increase with solute content (5psc = 
3400 nfim). Good predictions were obtained with this experimental value. 
As for the solute diffusion coefficients and the interface free energies, they 
were directly deduced [11] from an atomic model previously developed for 
Al-Zr-Sc alloys [15,18]. The good agreement obtained between our resistivity 
simulations and experimental data allows then us to stress the correctness of 
the AI3SC interface free energies used in our simulations. For temperatures 
ranging between 200 and 500° C, this interface free energy corresponding to 
coherent AI3SC precipitates is varying between 119 and 105 mJ.m~^. 
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One of the key assumptions of our simulations is that all clusters contribute 
to electrical resistivity and that each cluster contribution is proportional to 
its area. Although this assumption may look crude, it leads to quantitative 
predictions. In particular, it manages to reproduce the resistivity excess and 
its slow decrease during coarsening. This excess resistivity mainly arises from 
large clusters contributions whereas the solid solution contribution can be ne- 
glected. As a consequence, resistivity measurements during coarsening do not 
really allow to follow the solid solution concentration. In particular, resistivity 
do not obey LSW theory in a phase-separating system like supersaturated Al- 
Sc alloys at variance with the solid solution concentration. This involves that 
one cannot deduce from these measurements correct values of the precipitate 
interface free energy or of the solute diffusion coefficients. On the other hand, 
solubility limits obtained from resistivity experiments are correct as both the 
normalized resistivity and the solid solution concentration tend to the same 
value for long enough annealing times. 



A Reinterpretation of R0yset experimental data 



Usually, one uses electrical resistivity measurements to define the fraction of 
precipitated solute 

Pit) - P° 

where p{t) is the resistivity of the solid solution measured at time t and p° and 
p'^^ are respectively the initial and equilibrium resistivities. If the resistivity 
was truly proportional to the solute concentration like it can be assumed in 
an under-saturated solid solution, this definition would be equivalent to the 
one obtained from considering the solid solution concentration instead of the 
resistivity. But in a phase-separating system like supersaturated Al-Sc solid 
solution, both definitions lead to different quantities. 

In their work [10], R0yset and Ryum used Eq. lA.ll to define the precipitated 
fraction. Assuming that resistivity varies linearly with in the coarsening 
stage, extrapolation of data to infinite time was used to estimate the equilib- 
rium resistivity p'^'^. This leads to a correct value of p'^'^ for the higher annealing 
temperatures but not for the lower ones (T < 330°C). Indeed, as we previously 
saw, one cannot assume that resistivity follows the same time dependence as 
the solid solution concentration due to the large cluster contributions to elec- 
trical resistivity. Therefore p'^'^ cannot be obtained from such an extrapolation 
to infinite time, and it has to be calculated from the equilibrium Sc solubility in 
aluminum and the linear relation [6,21] between resistivity and concentration 
existing in dilute solid solutions leading to 

= p° 1 + 6ps.xll (A.2) 
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Fig. A.l. Evolution with time of the precipitated fraction Xp in an aluminum solid 
solution of composition Xg^, = 0.12 at.%. Xp has been deduced from experimen- 
tal resistivity measurements [10] considering the equilibrium resistivity is reached 
(R0yset data) or not (after correction) for the longest annealing time. 

The precipitated fraction Xp obtained calculating p'^^ in this way differs from 
the one obtained by R0yset and Ryum who assumed that p'^'^ corresponds to 
the resistivity measurements for the longest annealing time (Fig. lA.ll) . This is 
a clear manifestation of the large cluster contributions to electrical resistivity. 
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